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Abstract 

The Einstein's field equations of FRW universes filled with a dissipative 
fluid described by full theory of causal transport equations are analyzed. 
New exact solutions are found using a non-local transformations on the 
nonlinear differential equation for the Hubble factor. The stability of the 
de Sitter and asymptotically friedmannian solutions are analyzed using 
Lyapunov function method. 
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1 Introduction 



Quantum effects played a important role in the early Universe. For instance, 
vacuum polarisation terms and particle production arise from a quantum descrip- 
tion of matter. Some other processes capable of producing important dissipative 
stresses include interactions between matter and radiation [||] and different com- 
ponents of dark matter 0. Phenomenologically, these processes may be modeled 
in terms of a classical bulk viscosity 

A relativistic second-order theory of non-equilibrium thermodynamics was 
developed in Q. This formulation generalizes the expression of the entropy flux 
by the inclusion of quadratic terms in dissipative non-equilibrium magnitudes. 
Another formulation of this theory, called Extended Irreversible Thermodynamics 
has been made in [H] M. There, non-equilibrium magnitudes are introduced and 
a generalized Gibbs equation including the dissipative magnitudes is employed. 

Recently, qualitative analysis of homogeneous cosmological models using the 
full causal transport theory derived from the extended irreversible thermody- 
namic @, H has been made to investigate the stability of the particular solu- 
tions of the Einstein's field equations. Several authors studied the stability of non 
flat causal cosmological models using non usual equations of state p/9 2 = Pox 1 , 
C/Q — (ox™ 0- They are able to get an autonomous system of differential equa- 
tions in a phase space (x, y) where x = 3p/0 2 and y = 911/ 9 2 . Other authors |§ 
have studied flat models with the widely used equations of state p — (7 — l)p 
and ( = ap m . They have got an autonomous system in the phase space (H, H). 
In both studies, always the linearization technique near stationary solutions in 
the phase space have been used. Some other papers where the full version of the 
transport equation has been used are || |10| [JTlJ] [12[ . 



In this work, together with another one new exact cosmological solutions 
of full causal models are found and the asymptotic stability of solutions present in 
non-causal and truncated causal models is investigated using the direct method of 
Lyapunov. Few authors have used this method to analyze the stability of causal 
models though it provides useful information of the dynamical behaviour of the 
system, not only near the stationary solutions but also far away from them |13| 



H| ||1 6[ . In the present paper we assume that the cosmic fluid is described locally 



by the equation of state p = nT, while in the paper |]13[ a power-law relationship 
between the equilibrium temperature and energy density was assumed. 

This work summarizes as follows. In section II we state the basic equations 
of the model, they are solved in section III for m = 1/2, the two-parameter 
families of solutions are classified and their asymptotic behaviors are analyzed. 
Also, the thermodynamics of these solutions is considered. In section IV we study 
the stability conditions for de Sitter and asymptotically Friedman solutions for 
1/2. Finally in section V we summarize the main conclusions of the paper. 
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2 Hydrodynamic description for the cosmic fluid. 
Basic equations 



For the fluid approach to remain valid, the mean interaction time t c of the fluid 
particles must be much less than the expansion rate. This requires, t c <C H^ 1 , 
where H is the Hubble expansion rate. If the fluid is in or close to equilibrium, 
then this restriction ensures that the fluid has a well-defined local temperature. 
For a cosmic fluid without shear viscosity or heat flow, the energy momentum 
tensor is 



Tab = (p + p + U)u a u b + (p + Il)g ab , (1) 

where u a is the four velocity, p the energy density, p the equilibrium pressure and 
II the bulk viscous pressure. In an expanding fluid, the dissipation due to II leads 
to a decrease in kinetic energy and therefore in pressure, so II < 0. The particle 
flow vector N a is given by N a = nu a where n is the particle number density. In 
absence of vector and tensor dissipation, the entropy four-flux will be of the form 
S a = S e ffn a where S e ff is the effective non-equilibrium specific entropy which 
must be positive. 

In the Israel/Stewart theory 

S eff = S T —Ii 2 (2) 

n 2nT( v ; 

where r is the relaxation time, T is the local-equilibrium temperature and ( is 
the coefficient of bulk viscosity. The conservation laws N a a = and T a \ = 
imply 



n + 3Hn = and p + 3H(p + p + IT) = 



(3) 



respectively, where 3H = u a . a is the fluid expansion and h = n )a u a . Combining 
equations (§) with the Gibbs equation for the local-equilibrium variables S and 
T 



TdS = d 



nj \nj 



we get 



S 



3HU 

nT 



From (Q) and ([3p we find 

na 

D .„ — 



C 2 I (T 



(4) 
(5) 

(6) 
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for the entropy production rate. The simplest way to guarantee the second law 
of thermodynamics S a a > is establishing the transport equation 

n + rri = -3(H- -rll [3H + ~ - - - . (7) 

for II, leading to S? a = II 2 /QT. 

We need to specify equations of state p = p(T, n) and p = p(T, n) for the 
cosmological fluid. In order to have a simple model, we assume that the cosmic 
viscous fluid may be locally described by an equilibrium equation of state that 
has the approximated form of an ideal gas and the 7-law: 



p = nT 

p = -J— nT. (8) 
7 - 1 

Note that both T and n are independent variables and only in the equilibrium 
limit the particle number density depends exclusively on T. 
From (|) and © the expression for the temperature is 

T = *^1r 3 p. (9) 

We specialise to a flat FRW universe, where the field equations are 



H 2 = ^p, 



1 

3' 



3(H + H 2 ) = + 3p + 3IT), (10) 

where H = R/R is the Hubble factor, R(t) the cosmic scale factor of the 
Robert son- Walker metric, 87rG = c = 1 and an overdot denotes differentiation 



respect to time t. The relaxation time r is related to the bulk viscosity by [14 



r = -T- ( n ) 

where v is the dissipative contribution to the speed of sound. By causality v 2 < 
2 — 7. The bulk viscosity coefficient is often taken to be C = ap m where a is a 
positive constant. 

An other interesting restriction is that, since a dissipative expansion is non- 
thermalising, the rate for some interaction (~ t _1 ), which is crucial to main- 
taining equilibrium, must remain lower that the expansion rate H, so that the 
fluid cannot adjust sufficiently rapidly to the changes in temperature induced by 
expansion 
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Finally, using (|7|), (|J), ( [XQ|) and (|TTD we get, for the evolution of H, the 
equation 



HH + flHH 3 ~ 2m - 2H 2 + ^Jll ( 7 H i-2m _ 3 
d 2 \o 







(12) 



where 5 = a3 m 1 , and we are assuming that the Universe is expanding, that is 
H > 0. Using (|), (|) and (|I0|) we get, after integrating 



S(t) = So + 



7-1 



hi(H 2 R 31 ). 



(13) 



3 Exact solution for m = 1/2 



The Einstein's field equations of FRW universes filled with a dissipative fluid 
described by both the truncated and non-truncated causal transport equations 
can be linearized and solved for a power-law fluid when m = 1/2 |13| Hl5| . In 
what follows, we shall find the general solution for a fluid that obeys (H). In this 



case setting m = 1/2 in (|12|) we find 



v 2/ y ■ o 3aH 4 ( 7 
HH + —-HH 2 -2H 2 + 



2 a- 3 



0. 



(14) 



With the change of variables H = 1/s, t = dcr/v 2 ^, equation ([14]) becomes into 

d 2 s 1 ds fi 



da 2 s da s 

where /1 = (36 /2v 2 7 )(3<5 - 7). 

Equation fllED is a particular case of the most general equation 



(15) 



s + «i/(s)s + a 2 f(s) J f(s)ds + a s f(s) = 0. 



(16) 



which linearizes through the nonlocal transformation [|15| [|17]] [|T^] |]T9| 

z = J f(s)ds V = J f(s)dt (17) 

becoming 

z" + ol\z' + a%z + «3 = (18) 

where the ' indicates derivative with respect to the new variable 77. In our case 
a.\ — 1, CH2 = 0, 03 = /i, f(s) = 1/s and the nonlocal transformation is then 



z = Ins 



w 2 7 f dt 
T 



(19) 



Thus, integrating 



z{rj) = ci — fir] - c 2 e v (20) 
we obtain the general solution of ([14]) in parametric form 

H(rj) = H exp {jnq + c 2 e-") . (21) 

At = — ^— / d V exp (-/^ - c 2 e~ r ') . (22) 

where At = £ — to, and Hq, c 2 , to are arbitrary integration constants. The scale 
factor can also be found in terms of 77 as 

J R(r ? ) = J R exp(4-r/) (23) 
\v ^7 / 

where i?o is another integration constant. In this case ( ~ H and r ~ 

Below, we show explicitly exact solutions and the asymptotic behaviour of 
the general solution. 

3.1 General solution for fi = 1 



In the particular case fj, — 1, Equation (|T6|) has an additional form invariance 



symmetry. Then, introducing the transformation [T8 



- 1 <5 u> 1 

S = = 'v^Jw^dt ( ^ 

into (^) we find w"(a) = and w;(cr) = Cicr + c 2 . Finally, the two-parameter 
family of solutions is given by 

R(t) = R (c ± In I At I )- 5 /'" 2 \ (25) 

where 5 = yy + a/7 2 + 8f 2 7^ /6, and c is an arbitrary integration constant. There 
exist solutions that expand so quickly that the scalar curvature diverges when the 
scale factor diverges and they become infinite in a finite amount of proper time. 
Some of these singular solutions, those with minus, have a finite time span and 
particle horizons. Also, we have the time reversal of the above solutions. These 
with plus lead to a shrinking universe. In both cases, the dynamical information 
is enough in order to select the optimal physical solutions. Non-thermalising 



condition requires in this case v < ^^l + 2/3j 
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3.2 General solution for /i = 



We note first that in this case (7/3 = 5) equation ( |12|) admits a de Sitter stable 
solution H = Ho, where Hq is arbitrary ||. On the other hand, for H 7^ Hq, 



equation (|TJ) can be integrated twice by quadratures 



At 



dH 
H 2 



In 



H 



3v 2 J H 2 \'" H , 
and we find these two-parameter families of solutions in implicit form 



At 



1 



At 



3v 2 H 
1 



Ei(ln^ 



3v 2 H 



Ei 



Ro, 



for H/H > 1, and 



At 



At 



1 n ./ 
— - — Ei - In — 

3v 2 H V H 



3v 2 H 



Ei 



R_ 

Ro 



(26) 

(27) 
(28) 

(29) 
(30) 



for H/Hq < 1, and where Ei and Ei are exponential integrals. Solutions of 
both families have an asymptotically de Sitter stage for large times, while their 
respective initial behavior (singular or nonsingular) follows those of the solutions 
with n = 1. Non-thermalising condition requires here v < l/v3 and the second 
law of thermodynamics does not impose additional restrictions. 



3.3 Particular solution for C2 = 

If c 2 = in Q2g) then H(rj = 0) = H and R{rj = 0) = R and we find that 
H = 2/ (3(7 — 35) At) and the scale factor is given by 

R(t) = RqA^ (31) 

where (3 = 2/3 (7 — 35). This solution is also found if we set H — (3/ At in (|TJ). 
If (3 > 0, that is, for 5 < 7/3, the Universe starts from a Big-Bang singularity 
and then expands, while for 5 > 7/3 the scale factor diverges in a finite amount 
of proper time. For 5 = 7/3, we have \i = and z = c\. In this case R = Roe Hot 
is the de Sitter solution. Non-thermalising condition and S e ff > 0, at late times, 
require that for v < l/\/3, a must be in the range v / 3'W 2 7 < ce < 7/V0. 
For 

/ 1 [1 2~ ,7 2^3 7 

v < — = Wl and — = < a < — = 

v^V 3 7 v/3 9 y/3 
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there exist inflationary thermodynamically consistent power-law solutions given 
by ©. For 

R 2 ^ ^ 7 2x/3 

V of 7 < a < — = — 

V3 9 

non-inflationary friedmannian thermodynamically consistent solutions given by 
( |31~D exist. We conclude that whenever 

1 /; 2" 

t> < — =4/1 

inflationary power-law and non-inflationary scenarios are thermodynamically 
consistent. 

For (3 = 1/2, that is a = V3(3j - 4)/9 for 4/3 < 7 < 2, the universe behave 
as if were filled by a perfect fluid and were radiation dominated. For (3 = 2/3, 
that is, a = ^3(7 - l)/3, the universe behave as if were filled by a perfect fluid 
and were matter dominated. 



3.4 Analysis of the solutions 



Here we analyze the behavior of the general solution (21)- (p2|). We begin by 
studying its asymptotic behavior near the singularities and for large times. 

For H — > 00, s — > and z — > —00. This occurs in two cases. If fi > 0, z ~ —fir] 
for 77 — > +00. Then, using ([19]), we find that H ~ e M?? = — l/(/icr). So, there is 
a two-parameter family of solutions with an explosive singularity at a finite time 
whose leading behavior is given by (|3TD . Also, H diverges when 77 — > —00 and 
C2 > 0. In this case, H ~ exp (c2e~ v ) and there is a two-parameter family of 
solutions representing a Universe beginning to from a Big-Bang singularity as 

R(t) ~ I In |A/;|r 31 ' 2 for At -> (32) 

The opposite extreme is H — > 0. In this case s — > 00 and 2 — > +00. If 
H < we obtain 2 ~| /i | r\ for r7 — * +00. Using the ([19]) again, we obtain 
H = e"!^ = 1/(| /i I a). So there exists a biparametric family of asymptotically 
Friedmann solutions with leading behavior given by (|31|). Furthermore, these 



solutions are stable, and fl3ll) is an attractor. 

When i] — > —00, two different behaviors occur depending on whether C2 < or 
C2 = 0. In the first case, H ~ exp (c2e~ v ) and there is a two-parameter family of 
solutions representing a Universe beginning to evolve as 

R(t) ~ I 1b |At|r 3 ^ for t -> -00 (33) 

In the second case the solution is (J5T]) . 

The two-parameter families of solutions can be classified according to their 
number of singularities. 
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There is a family of singularity- free solutions which begin at t = — oo and 
their behavior in the far future is Friedmann or have a final de Sitter stage. 
These solutions do no have particle horizons. 

There is a family of solutions that begin as in the previous case but end at an 
explosive singularity. These solutions do not have particle horizons. Also there 
is a family of solutions that begin at a Big-Bang singularity and have a final 
Friedmann or de Sitter behavior. 

Finally, there is a family of solutions with a finite time-span that begin at a 
Big-Bang singularity and end at an explosive singularity. 

3.5 Thermal consequences 

Inserting the general solution (|2lD, (p2|) in (§) 



T = 3- — -R 3 H$exp 
n 



^- (1 + 35 - 7) 77 + 2c 2 e~ 



(34) 

v ^7 

we can investigate the evolution of the temperature for the different families 
of solutions studying the limits near the singularities and for large times. For 
t — ► oo, the temperature vanishes provided that a < v^3(7 — l)/3, so the scale 
factor grows asymptotically slower than t 2 ^ 3 . The above condition implies that 
fx < — (7 — 1) / (2^v) < 0. This shows that, although a power-law inflationary 
scenario is possible, the temperature grows unboundly in the far future. For 
evolutions with a singularity in the past, we have a hot Big-Bang beginning of 
the Universe because the temperature diverges at the singularity. On the other 
hand, for evolutions ending at an explosive singularity, the temperature diverges 
as a negative power of the proper time interval to the singularity. There are two 
two-parameter families of solutions which have a cool beginning with R = in 
the far past. These contain a set of solutions that have a cool ending provided 
a < v / 3(7 — l)/3, otherwise they have a hot final stage. With a fine-tuning of 
the parameters which corresponds to solutions with the asymptotical behavior 
R ~ t 2 / 3 when t — > 00, T — > 3RIHq (7 — 1) /n . In this particular case, the 
solution (|31"D has a constant temperature. 



4 Stability analysis for m^l/2 



We will make use of the method of the Lyapunov function [|16j to investigate the 
asymptotic stability of the de Sitter and asymptotically Friedmann solutions that 
occur in the noncausal and truncated causal models [15[], as well as in the full 
causal model a with power-law equation of state [T3 . 
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4.1 Stability of the de Sitter solution 

Equation fll~2]) admits a de Sitter solution [§] 



/or\ l/(l-2m) 

tf = f (m ^ 1/2). (35) 

To study its stability we make the change of variable if = s _1 in ( |i~2"D which 
becomes into 

S + ^s 2 - 2 s - ^^ 2m - 2 + ^ = (36) 
multiplying by s one may rewrite (^) as 

where the term s 2 /2 represents the kinetic energy term of the system and V(s) 
the potential energy term. The dissipative term D(s, s) is negative definite for 
all s. One may calculate these terms and finds 

K (S = ; 1 Ins , 

v ' 26 2m - 1 2 

jD ( s ,s) = -4 I sV m - 2 . (38) 
o 

The potential V(s) has a unique minimum for s > at so = -f^o" 1 for m < 1/2, 
while for m > 1/2 it becomes a maximum. Thus we find that an exponential 
inflationary regime is asymptotically stable for t — > oo for any initial condition 
H > provided that m < 1/2 (node or stable focus), but this regime becomes 
unstable for m > 1/2 (saddle point). This result improves over previous studies 
based in small perturbations about the solution 



4.2 Stability of the asymptotically Friedmann solution 



For m > 1/2 it is easy to check that (|T2| ) admits a solution whose leading term 
is 2/(37t) for t — > oo. This suggests to make the change of variables 

H = — *— , t 2m - x = z (39) 



in (|1^) which takes the form 

■ D(u,u',z) (40) 



^u' 2 + U(u,z) 



dz 

where ' = d/dz. We consider that u lies in a neighbourhood of Uq = 3 / ~y/2. Thus 
when z — > oo 
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,„ v l,2 7 /f 2 ™ Slii 2 '"- 1 ^ 1/1. 



for 



for m — 0, and 



B <"'"'^=-j(2^i)" 2mM ' 2+o G)- (43) 

As U(u,z) has a unique minimum at u for any m 7^ 1/2, and D(u,u',z) 
is negative definite for m > 1/2, we find that solutions with leading Friedmann 
behaviour R ~ £ 2 /( 3 t) when t — > 00 are asymptotically stable for m > 1/2, but 
become unstable for m < 1/2. This behavior has already been found in the 
truncated theory and the power-law model pi] 



5 Conclusions 

A cosmological model for a causal viscous fluid with an equation of state like that 
of an ideal gas is analyzed carefully. Our main result is that many qualitative 
features of this model are similar to those of the full causal model with a power- 
law equation of state for the temperature and the truncated model. 

Applying a non-local integral transformation to the nonlinear differential 
equation for the expansion rate it linearizes, and we find the general solution 
of this model for m = 1/2. Due to viscous effects, these new exact solutions 
contain a two-parameter family that avoid the singularity and do not have parti- 
cle horizons. This interesting evolution has already been found in the full causal 
theory, with a power-law equation of state for the equilibrium temperature JT3 



Among the set of solutions for some values of the parameters we have found 
explicit exact ones that give the qualitatively behavior of parametric solutions 
for \i > 0. Also, there are solutions that begin at a Big-Bang singularity and may 
represent a power-law inflationary stage or a matter or radiation dominated era. 

For m = 1/2, the splitting of the large time asymptotic behavior of solutions 
in terms of sgn(7 — 35) follows closely that of the power-law model and is the 
same as in the truncated model, which in its turn resembles the classification 
for the noncausal solutions. As in the power-law case, no family of solutions 
with extrema is found, and this appears to be the main consequence of the full 
theory. On the other hand, a logarithmic leading behavior appears near the Big- 
Bang singularity and in the remote past. This kind of evolution has not been 
observed previously for these models. We note that the rate of expansion in this 
logarithmic Big-Bang is larger than in a regular power-law one. 
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The asymptotic behavior of solutions for this model was analyzed using the 
method of Lyapunov for m 7^ 1/2. We have shown that friedmannian solutions 
are an attractor set of the biparametric family of solutions for m > 1/2, while 
they are unstable for m < 1/2. This information is not available via the stability 
analysis by linearization ||. So we have improved this analysis, obtaining more 
information away from the stationary solutions. Also we have found that an 
exponential inflationary stage is asymptotically stable for m < 1/2 while it is 
unstable for m > 1/2. Here also, the same structural behavior as in the power- 
law, truncated |2(| and non-causal models is observed. 

It has been conjectured that due to the use of causal truncated or non- 
causal theories, viscosity-driven exponential inflationary solutions are spurious 
0. There, a cosmological model containing a dissipative Boltzmann gas is stud- 
ied. We show here that this conclusion does not generalize to models described 
by another equation of state. 



Acknowledgments 

We thank to Diego Pavon for his interesting comments about the results 
of this paper. One of us (LPCh) wants to acknowledge the hospitality of the 
Departament de Fisica, Facultat de Ciencies, of the Universitat Autonoma de 
Barcelona, where part of this work was done. This was partially supported 
by project 3091/92 CONICET and by the Direction General de Investigation 
Cientifica y Tecnica of the Spanish Ministry of Education and Science under 
grant PB94-0718. 

References 

[1] Weinberg S 1971 Astrophys. J. 168 175 

[2] Pavon D and Zimdahl W 1993 Phys. Lett. A 179 261 

[3] Hu B. L. 1982 Phys. Lett. A 90 375 

[4] Israel W and Steward J M 1976 Phys. Lett. 58A 213 

[5] Pavon D, Jou D and Casas- Vazquez J 1982 Ann. Lnst. H. Poincare A 36 79 

[6] Jou D, Casas- Vazquez J and Lebon G 1993 Extended Irreversible Thermo- 
dynamics (Springer Verlag, Berlin) 



12 



[7] Coley A. A., van der Hoogen R J , Maartens R. 1996 Phys. Rev. D 54 1393 

Mendez V. and Triginer J. 1996 J. Math. Phys. 37, 2906 (1996) 

[9] Hiscock W A and Salmonson J 1991 Phys. Rev. D 43 3249 

[10] Zakari M and Jou D 1993 Phys. Rev. D 48 1597 

[11] Romano V and Pavon D 1994 Phys. Rev. D 50 2575 

[12] Maartens R 1995 Class. Quantum Grav. 12 1455 

[13] Chimento L P and Jakubi A S 1996 Class. Quantum Grav. (in press) 

[14] R. Maartens, in Proceedings of Hanno Rund Workshop, to appear; |astro- 



ph/9609119 



[15] Chimento L P and Jakubi A S 1993 Class. Quantum Grav.10 2047 

[16] Krasovskii N N 1963 Stability of Motion (Stanford University Press: Stan- 
ford) 

[17] Chimento L P and Jakubi A S 1996 Phys. Lett. A 212 320 
[18] Chimento L P 1996 J. Math. Phys. (in press) 

[19] Aguirregabiria J M and Chimento L P 1996 Class. Quantum Grav. 13 3197 
[20] Pavon D, Bafaluy J and Jou D 1991 Class. Quantum Grav. 8 347 



13 



